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We study analytically as well as numerically the dynamics of a quantum map near 

■ a quantum resonance of an order q. The map is embedded into a continuous unitary 
transformation generated by a time-independent quasi-Hamiltonian. Such a Hamiltonian 
generates at the very point of the resonance a local gauge transformation described the 
unitary unimodular group SU(q). The resonant energy growth of is attributed to the zero 
Liouville eigenmodes of the generator in the adjoint representation of the group while the 
non-zero modes yield saturating with time contribution. In a vicinity of a given resonance, 
the quasi-Hamiltonian is then found in the form of power expansion with respect to the 
detuning from the resonance. The problem is related in this way to the motion along a 
circle in a (q 2 — l)-component inhomogeneous "magnetic" field of a quantum particle with 
q intrinsic degrees of freedom described by the SU (q) group. This motion is in parallel 
with the classical phase oscillations near a non-linear resonance. The most important role 
is played by the resonances with the orders much smaller than the typical localization 
length, q <C I. Such resonances master for exponentially long though finite times the 
motion in some domains around them. Explicit analytical solution is possible for a few 
lowest and strongest resonances. 
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1 Introduction 



Classical canonical two-dimensional maps originated from the Poincare sections in the phase 
space has played an exceptional role in the establishing of our understanding of the origin 
and properties of the dynamical chaos [I], 0- Formally, they correspond to non-conservative 
Hamiltonian systems with one degree of freedom driven by instantaneous periodic kicks. The 
phase plane of such a map is, generally, very complex and consists of intimately entangled 
domains filled by regular and chaotic trajectories. The chaotic domains remain isolated from 
each other if the driven force is weak, but they join and global chaos appears when the strength 
of the force exceeds some critical value . After that, regular motion survives only inside isolated 
islands of the phase plane, where phase oscillations near the points of mainly low order non- 
linear resonances take place, whose areas diminish with the strength growing. 

Existence of the chaotic domains signifies absence of a global analytical integral of the 
motion. At the same time, there exist in the regions of the regular motion approximate local 
integrals which can be chosen in many different ways. It should seem, however, that the most 
convenient and physically sounding choice is that of the quasienergy integral which is directly 
linked to the periodicity of the driving force. In such an approach ||, the regular regime near 
a non-linear resonance is juxtaposed with the continuous evolution described by a conservative 
effective quasi-Hamiltonian function. Corresponding canonical Hamilton equations generate 
continuous trajectories on which all phase points of the original map lie. The quasi-Hamiltonian 
is found by perturbation expansion with respect to a small parameter which can, in particular, 
be the closeness to the resonance. 

After the quantum extension of the canonical maps had been suggested in M , the quantum 
maps were widely used as informative models of quantum chaos. Amongst them Chirikov's 
standard map, i.e. the periodically driven planar rotor, proved to be the most economic, 
fruitful and popular. The unitary Floquet transformation U which evolves the QKR wave 
function ip{9) over each kick period is given by: 



U = U r ■ Uk = exp ^— -Trh 2 J ■ exp (—ikcos8) (1) 

and consists of successive kick transformation Uk with the strength k and a free rotation U r 
during the time T. Here m = —id/d9 and we put H = 1. The standard map provides a 
local description for a large class of dynamical systems. In particular, there exists a tight and 
remarkable analogy || between discovered in |4| dynamical suppression of chaos in QKR and 
Anderson localization in quasi ID disordered wires. The diffusive growth of the QKR energy 
turns out to be restricted to a certain maximal value because of dynamical localization in the 
angular momentum space. 

At the same time, there exit some important features of the QKR dynamics, namely so 
called quantum resonances, which have no counterparts in the disordered systems (see || and 
discussion in [0, H). At a fixed value of the kick parameter k, special resonant regimes of motion 
appear || for everywhere dense set of the rational values = T/4tt — p/q where the integers 
p and q are mutually prime. Under these conditions the rotator regularly accumulates energy 
which grows quadratically in the time asymptotics fP2| , [TjJ. Both the restricted diffusion and 



the quantum resonances were experimentally observed in the atom optics imitation of the QKR 
reported in |J . Quite recently, the regime of quantum resonances re-appeared in a new aspect 
in connection with the electron scattering with excitation of the Wannier-Stark resonances 
PI PH. 
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Our results have already drown much interest of the authors (private communication). 

The interesting and important problem of the impact of the quantum resonances on the QKR 
dynamics and the interplay between resonant and diffusive regimes is still far from a satisfactory 
understanding. Investigation of this problem is the main goal of the present paper. We develop 
a general approach to the problem of the motion in a vicinity of a quantum resonance with 
an arbitrary order q. Generally, the influence of a quantum resonance depends on the relation 
between the order q and the localization length I in the angular momentum space. We show that 
the most important role is played by the resonances with q <C I- for finite regions around them 
the quantum motion is explicitly shown to be regular and dominates the motion for all values <; 
inside these regions - the resonance widths. More precisely, the motion is well described, during 
large though finite times, by a time-independent effective quasi-Hamiltonian with one rotational 
degree of freedom and with a discrete spectrum. Such a motion is in parallel with the classical 
phase oscillations near a non-linear resonance. On the contrary, when the resonance order is 
large enough, q ^> / the resonant quadratic growth appears only in the remote time asymptotic 
and for lesser times the motion reveals universal features characteristic of the localized quantum 
chaos. 

In sec. 2 we explain the concept of the effective quasi-Hamiltonian on which our approach is 
based. The power expansion of the quasi-Hamiltonian near a quantum resonance is constructed 
in sec. 3. In this section we also consider analytically and numerically two strongest boundary 
quantum resonances with q = 1, 2 and their classical limits. Two more strong resonances are 
investigated in the next sec. 4. Contrary to the boundary resonances, they disappear in the 
formal limit h — > 0. General consideration of a resonance of an arbitrary order is presented in 
sec. 5. At last, the problem of convergence of our expansion is discussed in sec. 6. 

2 Quasi-Hamiltonian 

Evolution of the QKR wave function for n kicks is given by n successive repetitions U (n) = U n 
of the Floquet transformation ([!]). Being unitary, the latter can be expressed in terms of a 
hermitian operator 7i as U = exp(— iH). Let us now consider continuous unitary transformation 
U{t) = exp(-iHt). According to such a definition, the wave function ip(9;t) = U(t)ip(9;0) at 
the integer moments t = n coincides with the quantum state of the map (0). On the other 
hand, the function ijj(8]t) satisfies standard Schrodinger equation with the time-independent 
Hamiltonian 7i. Obviously, the very possibility of the formal construction described is based on 
the periodicity of the map. That is why we refer below to this operator as the quasi-Hamiltonian. 

Let |e) be the eigenvector of the Floquet operator ([[]), which belongs to an eigenvalue e~ le . 
Then the quasi-Hamiltonian can be expressed as 



where the sum runs over the quasienergy spectrum {e} of the rotator. As usual, each quasienergy 
is defined up to a term multiple 2tt which results in corresponding ambiguity of the quasi- 
Hamiltonian (0). However, each time one can fix the operator 7i in the way most convenient 
for calculation. The ambiguity does not influence the evolution operator U(t) at integer mo- 
ments. As a rule, to get rid of the ambiguity we suggest continuity of the quasi-Hamiltonian 
with respect to parameters under consideration. In the coordinate representation, the quasi- 
Hamiltonian is an operator function of the pair of canonically conjugate observables 9 and m. 




(2) 



e 
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Since the operators m 2 and cos 9 do not belong to any finite algebra, the operator TC cannot 
generally be found in a closed form. Rather, it is expressed as an infinite sum of successive com- 
mutators. More than that, one anticipates extremely non- uniform dependence of the operator 
TC on the parameters <j and k. 

However, the problem simplifies enormously if the condition of a quantum resonance is 
fulfilled. As it has been shown by F. Izrailev and D. Shepelyansky at the point of a 
quantum resonances with the order q the Floquet operator can, generally, be presented as 
a q x q matrix. This matrix turns out to be (see below) a local gauge transformation from 
the unitary SU(q) group generated by a hermitian matrix TC^ res \9) which depends only on 
the angle and does not contain the angular momentum. Owing to this fact, the problem of 
calculation of the matrix 7v res ) becomes as simple (or complicated) as that of diagonalization 
of a g-dimensional unitary matrix. The latter can be carried out analytically if the matrix order 
does not exceed 4. 

Dependence of the quasi-Hamiltonian on the angular momentum recovers out of the points 
of quantum resonances. In some domain near a given quantum resonance with an order q 
this dependence can be found in the form of a power expansion over the detuning from the 
resonance. This expansion turns out to appear in the form of series, in particular, in powers 
of the angular momentum (see eq. (|105| ) below) with the resonance matrix TC^ res \9) being 
the zero-order term in the series. Actually, such an expansion is quite a formal one. The 
question of convergence by no means is trivial. At best, the series is of only asymptotic nature. 
Nevertheless, we shall see below that a few its first terms give surprisingly good description of 
the evolution during a very long time. 



3 The Boundary Resonances 

3.1 Regularity Domains and Quasi-Hamiltonians 

In the simplest case q — 1, the rotation operator U r is equivalent to the identity and the 
QKR evolution during a time t is described by the unitary phase transformation e~ lvt ; v(9) = 
k cos 9. This transformation parametrically depends on the angle 9 and therefore has continuous 
eigenvalue spectrum e(9) = v (9). By the moment t, an isotropic initial state, which we suggest 
throughout the paper, evolves into the wave function 

ip(9;t) — exp(— ikt cos 9) (3) 

V 2ir 

with ~ kt harmonics. The natural probe of the number of harmonics is the angular momentum 
operator rh = —id/d9 whose time evolution obeys the linear law 

m{t) =e ivt me- ivt = m + e ivt [m,e- ivt ]_ = m-v't-, v' = dv/d9. (4) 

This yields the quadratic growth of the kinetic energy E(t) of the rotor 

E(t) = l -{[m{t) - mf) = \{{v'f)t 2 = ^t 2 (5) 

with the resonant growth rate r = k 2 /A. Here and below the angular brackets denote averaging 
over the angle 9. 
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Let us now consider a vicinity k = T—Airp of a main resonance q = 1. The time-independent 
quasi-Hamiltonian in the vicinity is introduced by representing the Floquet operator in the 
form 



U Pi i(k) = exp y—-Km J exp(-ivt) = exp ^ Hi(k) j (6) 

with 

Hx(k) = kv + k 2 Q(k). (7) 
It follows that the operator Q(k) must satisfy the condition 

exp (^ l -Km^j = T*exp {-*«^ dr Q(k; -r)| , Q(k; -t) = e~ iVT Q(k) e ivr (8) 

where the symbol T* stands for the anti-chronological ordering. Suggesting the operator Q{k) 
to permit expanding in power series Q(k) = + k,Q^ + ... over the detuning k, one comes 
to successive relations 



Jo 



drQ^\-r) = -m\ (9) 



i l Tl 

| drgW(-r) = --J d n J dr 2 [Q^(-n), Q {0 \-r 2 )] _ (10) 



and so on, which allow to find the quasi-Hamiltonian (|7|) up to desired accuracy. Eq. (j9|) 
implies that the evolution described by the map ([!]) is smoothed in such a way that the overall 
effect of the lowest order correction k 2 Q^ within one period is identical to the kinetic energy 
operator /C = J 2 /2; kiti = J. Being a small factor in front of the angle derivative, the detuning 
k plays here the role of the dimensionless Planck's constant while J is the angular momentum 
operator in the units chosen. With only the two first corrections (P,|10D being retained the 
quasi-Hamiltonian acquires the form 

Hi = \jF 2 {9)J + I U + F (B) (11) 

of the Hamiltonian of a generalized pendulum. The periodic functions Fi(9) depend on the 
angle via the kick potential v(8). In the lowest approximation 

F 2 (6) = l; F 1 (6) = -^v>; F (6) = nv + ^(v')\ (12) 

when the next correction adds 

6F 2 (6) = -^v; 8F X {6) = ^/v 7 6F (6) = -^Wfv - ^v. (13) 

In further corrections higher powers of the operator J also arise. 

Uprising of the angular momentum operator J in (JTT|) drastically changes the eigenvalue 
problem. The angle 9 is a quantum-mechanical coordinate operator in this problem, and the 
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spectrum of the quasi-Hamiltonian Ji\ as well as of the Floquet operator becomes discrete be- 
cause of the periodic boundary condition. The main effect of the term quadratic in J consists 
in cutting off the unrestricted kinetic energy growth. This is well seen in Fig.l Already the 
lowest correction (dotted line) describes reasonably good the turn off from the quadratic res- 
onant growth, as well as the mean height of the saturation plateau. The next one (solid line) 
substantially improves the description and reproduces well also the details of quantum fluctu- 
ations in the plateau region. Influence of the further corrections, which contain, in particular, 
higher powers of the angular momentum J, remains weak in the finite domain Ak, the width of 
the resonance, where the angular momentum J in the plateau region is still sufficiently small. 

Before analyzing the conditions of validity of eq. ( |TT|) in more detail, we consider the reso- 
nance q = 2 because of a tangible similarity of the two cases. For q = 2, the rotation operator 

- 2 

U r = e~ l7rpm , where p is an odd number, has only two eigenvalues 1 and -1. Obviously, any 
periodic function ip+(9) (ip-(6)) which contains only even (odd) harmonics is an eigenfunction 
belonging to the eigenvalue 1 (-1). An arbitrary state ip(6) can be written down as the linear 
superposition ip = ip + + ip_ of these eigenfunctions. Therefore, in the Hilbert space of the pe- 
riodic functions the operator U r is isomorphic to the 2x2 Pauli matrix er 3 . One can actualize 
this isomorphism by representing the state ip(9) in the form of a two-component spinor 

m - . (i4) 

In this representation the differentiation operator rh which does not change harmonics' numbers 
grows into the diagonal matrix 

M -(o»)' (") 

while any periodic coordinate operator F{6) looks as 

F(0) = f + (9) I + a x (16) 

with I being the 2x2 unity matrix. At last, matrix elements of dynamical operators take the 
form 

02,1= j + d0#5(0)O(0,M)*i(0); j + dj9&(0)V(0) = l. (17) 

J — 7T J — TV 

The free rotation is described now by the matrix operator U r = a$ = e -1 ^^ 3-1 ) and the kick 
operator reads as = e~ lvai . Simple manipulations with Pauli matrices lead to the following 
expression for the resonant Floquet operator 

Ug s) =e^exp(-i| n-tr) (18) 

where the unit vector n = (0, sin v , cost" ). Up to the trivial phase factor, this is a spin-flip 
operator which belongs to the unitary unimodular group SU(2). Therefore, corresponding 
evolution in the continuous time t fully reduces to the spin rotation. 

In particular, evolution of the angular momentum is given by the equation 

AM(t) = exp (i| n- at) Mexp (—if n • <rt) — M = AM(t) • a. (19) 
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The vector AM(i) is easily calculated by making use of the formula 



exp ^— i— n • <j tj = cos(7rt/2) — in - cr sin(7rt/2). 

Simple transformations lead to the result 

AM(t) = -v [s sin(7rt/2) cos(?rf/2) + 1 sin 2 (7rt/2)] (20) 

where the unit vectors s and 1 are defined by 

n' = u'(0,cosu, -sin?;) = v's; 1 = [s x n] = (1, 0, 0). (21) 

The three unit vectors n, 1, s form an orthogonal basis in the 3-dimensional adjoint space of the 
577(2) group. The evolution ( p0|) is purely periodic in time, so that the kinetic energy 

E(t) = i([AM(t)] 2 ) = l -({v'f) sinV/2) = £ sinV/2) (22) 
does not grow but rather jumps between two values and A; 2 /4 when the time t runs over 



integer values [12 



The motion in a neighborhood of the considered resonance is described by the quasi- 
Hamiltonian matrix 



H 2 = /c|n ■ a + k 2 Q(k) (23) 

(compare with eq. (^)) where Q is a 2 x 2-matrix operator in the spinor space. This operator 
satisfies the condition (^) with rh substituted by the matrix M from eq. (IB). With the same 
accuracy as above, the quasi-Hamiltonian reads as in eq. (JTT|) again where the angular momen- 
tum J and the functions flT2"|) are replaced by 2 x 2 matrices. In the first approximation they 
are equal to 

F 2 (0) = /; Fi(0) = ~v' (l + |s) ■ o. V (ff) = ^n-cT + ^{v'f I. (24) 

The Hamiltonian (|TT|), ( p4]) describes the motion of a quantum particle with the spin 1/2 along a 
circle in an inhomogeneous magnetic field. The terms linear in the angular momentum J mimic 
a sort of the "spin-orbital" interaction. Fig. 2 (points) shows that quantitative agreement of 
this approximation with exact numerical simulations (circles) worsens rather fast. However, 
the second correction which can be represented in the compact form 

^ m - t, { J - { J - (i<"') 2 " - • CT L} + + lr + s<«'' 2 « (r - ') • " < 25 > 

noticeably improves the correspondence (crosses). The two branches correspond to even (starts 
from zero) and odd (starts from k 2 /A) kicks. 

At the first glance, an important difference is seen in the f-dependence of the functions Fi 
in eq. ([T2|) on the one hand side and in eq. (^4|) on the other. In the former case, the number of 
harmonics do not exceeds the power of the detuning k. This property holds also in the higher 
approximations. Likewise, the factors k and k are balanced in the similar way so that k always 
combines with k into the effective Chirikov's parameter K e = nk. Afterwards only positive 
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extra powers of k may remain. Therefore, the influence of the higher terms of the expansion 
can be expected to be weak when K e < 1. 

In contrast, the unit vectors n and s in eq . (|24j) in the most interesting case k ^> 1 contain 
~ k harmonics and their derivatives with respect to 9 are large. This leads to terms with extra 
powers of k in the higher corrections, which are not compensated by the small detuning k and 
enhance higher corrections. For example the term in (PB]), which is quadratic with respect to the 
operator J, contains uncompensated factor k. However, such terms turn out to be inefficient 
and cancel finally out due to the identity 



u: 



(res) 
P,2 



1 2 



exp 



-tir n • <x 



(26) 



so that the case q = 2 does not essentially differ from the main resonance as it concerns the 
role of higher corrections. This fact can be proven by disentangling the resonant part from the 
squared Floquet operator Up 2 (k) after which the non-diagonal cr-matrices disappear from the 
quasi-Hamiltonian. 

However, it is appreciably simpler to get rid of the non-trivial SU(2) algebra by separating 
evolution at only even and only odd kicks. It is then enough to smooth directly the squared 
Floquet transformation. Taking into account that in the Hilbert space of periodic functions 



cos 6* e 



cos#, one comes to the condition 



KM 



exp 



-—Km 
2 



exp 



-—K(m — v 
2 v ; 







= exp 


-2-H 2 (k) 

K 



The quasi-Hamiltonian is now found with the help of the Baker-Hausdorf expansion, 



(27) 



hie 



hi) 



5^- 



1 



24k 2 



fC — fC, [/c, /c] _ 



48k 3 



+ ... . 

(28) 



Here K = |J 2 and K = |(J — nv') 2 = |(J + Kksin9) 2 are the kinetic energy operators at 
the moments t = and t = 1 respectively. Since each commutator gives at least one power 
of the small detuning k, uncompensated factors k do not appear in the series. All four terms 
displayed in eq. ( p8"D are easily calculated explicitly though corresponding expressions are too 
lengthy for presenting them here. Fig. 3 demonstrates very good agreement of the evolution 
described by this quasi-Hamiltonian with exact numerical simulations. Only the even branch 
is shown. 

The q = 2 resonance of QKR is an example of the specific regimes of motion of periodically 
driven systems which are known as anti-resonances. The main feature of them is periodic exact 
recurrence (see eq. (^)) after a certain number of kicks. General consideration of the motion 
near anti-resonances is presented in [13]. The authors showed, in particular, that such a motion 
has regular classical limit. In fact, this is valid not only for q = 2 anti- resonance but also for 
the actual resonance q = 1 (see next subsection). More than that, we will demonstrate in 
sees. 4, 5 that domains of a regular quantum motion exist near all resonances with q <^ I. 
However, contrary to the two boundary resonances, the widths of other resonances diminish in 
the classical limit so that this motion has no direct counterpart in the classical standard map. 
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3.2 The Classical Limit 



In the both cases already considered the quantum fluctuations fade away when the parameter 
k increases. The quasi-Hamiltonians TCi and 7^2 appreciably simplify in the limit k ^> 1, k <C 
1; K e = nk = const and the functions Fi(Q) reduce to 

F 2 (9) = 1 - ^£ C os0; F 1 (9) = ^ sin - § sin 20, 
6 2 24 



K 2 \ K 2 K 3 

F (6) = K e (l M cos e cos 20 + cos 30 (29) 

1 240 / 24 240 V ; 



in the case q = 1 and to 

K K f *r>K 2 \ K 2 K 3 

F 2 (0) = 1 + — - — ( 1 + C os0 + ^cos20-^cos30, 

y ' 8 2 V 48 J 24 32 

/ 5K 2 \ K 2 K 3 K A 

FA 9) = — 1 + ^1 sin0 ^ sin 20 + ^ sin 30 - ^ sin 40, (30) 

w 2 V 48 / 8 96 96 ' v ; 

F = e - cos 20 ^ cos 40, 

v ; 8 192 

when q = 2. The Planck's constant disappear and corresponding quasi-Hamiltonians pass into 
classical Hamilton functions which depend on the only effective parameter K e . Treating K e as 
the classical Chirikov's parameter, these functions coincide with those obtained in ref. || and 
describe the phase oscillations near the non-linear resonances respectively of the first and second 
harmonics in the classical standard map. Fig. 4 illustrates on the example of the resonance q=2 
the quantum - classics correspondence. The solid line is obtained by averaging 1000 classical 
trajectories with J(t = 0) = over isotropic initial angular distribution. Crosses and circles 
show the results of numerical simulations of the exact quantum QKR map for two different 
values of the kick parameter k. The effective classical parameter is kept fixed, K e = 0.1. Being 
notably substantial for k = 10, the deviations due to quantum effects become inessential when 
k = 100. Agreement at large times are improved by taking into account terms of higher powers 
in J in the expansion of the frequency of the non-linear phase oscillations. 

The effective parameter K e differs from the product kT by the resonant part {4irp/q)k. In 
this connection, the suggestion made in ref. ]14j is worthy mentioning that in the regime of 



quantum diffusion the classical Chirikov's parameter is re-scaled as K e = 2&sin(T/2). This 
suggestion proved to be in good correspondence with numerical simulations. In the main order 
with respect to the small detuning k the re-scaled value K e is equivalent to our K e . However, 
it is not quite clear whether the whole sine make sense in the domains of the regular motion. 
The quantum corrections to different structures in the quasi-Hamiltonians have different forms 
neither of which can be identified with the terms of expansion of K e over the detuning k. 

We see that for both resonances considered the domains of regular motion is estimated by 
the same inequality K e < 1, which insures the regular phase oscillations near corresponding 
non-linear classical resonances. 

Note in conclusion of this section that near the resonances q = 1,2 the motion does not 
depend on the integer number p. This is a special manifestation of the following general 
properties of the standard quantum map. First of all, in virtue of the identity e ±2mm =^> 1 
the rotation operator U r and, consequently, the Floquet transformation (|1|) are periodic in 



9 



the parameter q with period 1. Therefore it is enough to restrict oneself to consideration 
of the interval < q < 1. In reality, only half of this interval exhausts all independent 



possibilities whilst the mean kinetic energy E(t) is calculated [TJJ. Indeed, presuming that 



q = 1/2 + 5q; \5q\ < 1/2, one can easily see that the transformation 8q — > —8q; 9 — > tx + 9 is 
equivalent to the complex conjugation of the operator U and therefore does not change E(t). 
Consequently, the problem investigated is symmetric in q with respect to the point q = 1/2 
and the lowest resonances q = 1, 2 correspond to the ends of the principal interval < q < 1/2. 
That is why we refer to these resonances as the boundary ones. 



4 Lowest Resonances with no Classical Limit 

Contrary to the two boundary resonances of the previous section, those of them which lie inside 
the principal domain < q < 1/2 have no well defined limits when k — ► oo, K e = const. We 
consider in this section the next two ones, q = 3,4; (p — 1). They provide typical though still 
exactly solvable examples of the motion near a quantum resonance. The resonance q = 3 is 
linked to the group 577(3) whereas the other one is still belongs to the same simplest group 
SU(2) as in the case q = 2. Indeed, it is easy to see that the rotation operator U r = e~ l ^ m is 
equivalent to the transformation e^ " 3-1 -* of the spinor (|14D [|I~7j| . Because of especially simple 



structure of this group we start with consideration of the resonance q = 4. 

Using well known properties of the Pauli matrices we easily find for the resonant Floquet 
transformation 

Again, we omit a trivial phase factor of no importance. The periodic functions in the exponent 
are now defined as 

( cos^A 1 
w(v ) = arccos — — , n(v ) = , — (smv , — smt>, — cost?). (32) 

V72/ Vl + sin 2 v 

The function w(v) is the most important new element in comparison to the q = 2 resonance, 
which yields a linearly increasing term in the angular momentum evolution 

v' 

AM(t) = -w'n t - V2 5- [s sm(wt) cos{wt) + 1 sm 2 (wt)} . (33) 

1 + sin v 

Here 

s= . = -^=(cosv , — cosf , 2 sinf). 1 = [s x nl = -^=(1, 1, 0). (34) 

Vl + sin 2 w V2 V2 J 

As before, the prime denotes differentiating with respect to the angle 9; the three unit vectors 
n, s and 1 are pairwise orthogonal. As long as the angle 9 remains fixed, the contribution of 
the spin rotation is periodic. However, on the last step averaging over the angle has to be done 
which leads to 

E(t) =rt 2 + X (t) (35) 

with 



1„ „ 2n k 2 r . 2 sin 2 ^ k 2 ( 4 f 1 

- {(wY) = — / d9sm 2 9 =— = — 1 / dz ^7— 

2 U ;/ An J_ n l + sin 2 w 4 ^ tt J 1 + sm 2 {kz) 
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and 

, , k 2 r ,„ . 2n sin 2 wt 2k 2 f 1 , r- sin 2 [w(kz)t] 

x (t) = — d9sm 2 9- ^— = / dzy/l - z 2 ^ — '-^ . 37 

W (l + sin 2 i;)2 vr J [l + sin 2 (^)] 2 

The function x(t) fluctuates with time slowly approaching the value 



k 2 f 1 , VT 



7T 



2 



2 



Xoo(fc) = — / <fc- . — 2 • (3* 



[l + sin 2 (A;; 



branches (s = 0, 1, 2, 



After extracting the constant part, the function naturally stratifies into four smooth 

X (+) (t) t = 4s, 
*H(t) t = 4s 
X {+) (t) t = As 

t = 4s + 3. 

Here 



r ,x, . k 2 f 77 , cos2wt i \ / \ k 2 r , sm2wt . , 

X (+) W = i -/ ^ , 2 ; x H (t ) = — / d« , 2 (40) 

where the function w(v) = ir/4 — arcsin(cos?;/"\/2) is 27r-periodic with respect to v and changes 
from zero at v = to the maximal value ir/2 when v = it. Stationary phase calculation gives 
for x (k) the asymtotics const /y/i. Fig.5 presents an example of the function x{t) taken at 
integer values of t. 

Simplification is possible in some limiting cases. It is easy to see that r(k) rs k A /16 |T2| and 
Xoo{k) ~ A; 2 /4 when the parameter k <C 1. In this limit w(kz) « 7r/4 + k 2 z 2 /2 so that 

= j [cos(£; 2 t/2) J (k 2 t/2) + sin(A; 2 t/2) J^kh/2)] 
X { '\t) = j [sin(£; 2 t/2) J (kH/2) - cos(k 2 t/2) J^kh/2)] . (41) 
The symbol J y (x) stands for the Bessel function. For small times k 2 t C 1 we find 

^'w-tO -£*")• ^"n*- (42) 

while 

X (+) (*)=X ( - ) (*)«^-0 (43) 

if the time is large. In the most interesting case k 3> 1 simplification is achieved by averaging the 
fast varying factors in the integrals over z in eqs. fl36|) , (]37D . This gives r(k) ~ (a/2 — l)A; 2 /4y / 2 
which is in good agreement with numerical data though somewhat differs from the value k 2 /12 
given in |~2|. At last, Xoo(k) ~ 3A; 2 /16v / 2 in this limit. One sees that the q = 4 resonance 
admits of very detailed analytical description. 
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Near the resonance calculation of the lowest correction gives 

F 2 = I; F x = ~v' {mn + V2(n\ + nf) [ws + (1 - n 3 w)l]} • cr; 

'1 



■ l\1 



F = kw n ■ cr + — (v ) 



~nl + (nl + nl) 2 (l-n 3 w) 



(44) 



(compare with eq. (p4"D). Computation of the next correction though making no principle 
problems, turns out to be rather tedious and leads to quite cumbersome expressions. The most 
important contribution comes from the correction to the matrix F 2 

6W 2 (9) = -1(1 - n?)(t/) V*i + \i 2 ) ■ cr (45) 
4 o 

with the vectors fi^ given by 

f x = (m(3 - 5n?)(l - n?), 1 + - 5n|, l + 2n\- 5nf) , 

f 2 = (rim 3 (14 - 15n?), -nm 3 (4 + 15n?), -(4 + 7n? - 30nf)) . 

We drop here corrections to other matrices F whose influence is negligibly weak. In Fig. 6 the 
evolution generated by the quasi-Hamiltonian (solid line) is compared to the simulation of the 
exact QKR quantum map which is shown by points. The detuning is chosen to be about 1/2 
of its critical value after which the regime of regular motion breaks into diffusion. The energy 
E(t) scales with the first power of the detuning k in this case. It is due to the form of the 
zero-order resonant interaction which contains spin and has no classical limit. The theory nicely 
reproduces all details of the evolution up to very large times. 

The width of the resonance is now much narrower than in the case of the boundary res- 
onances. Indeed, ^-derivatives of the functions w(v) and n(v) appear in higher corrections, 
which are large if k ^> 1. Validity of the expansion deteriorates because of such contributions. 
One can roughly estimate the width suggesting that the influence of the second order correction 
should be relatively weak inside the resonant domain. This gives Ak oc k~ 2 which agrees with 
our numerical data. Outside this interval the expansion transparently diverges. Contrary to 
the boundary resonances, the region of regular motion vanishes in the classical limit k — ► oo 
even if the condition K e = nk = const holds. 

To explore the regularity domain and adjacent area in more detail, we have fitted in Fig. 7a 
exact numerical data for the mean height E p i of the plateau as a function of the detuning k. 
Two qualitatively different regions are clearly seen: the regularity domain (k < 10~ 4 ), where 
the plateau is inversely proportional to the detuning, and the quantum chaos > 10~ d ) 

where the height is scattered around the generic value I 2 ~ k 4 . In the intermediate domain 
the higher corrections become increasingly important and the perturbation expansion fails. In 
Fig. 7b similar numerical results are displayed for the case q = 3. 

At the resonance point q = 3 the wave function is naturally split into three independent 
parts, t/>(0) = ipi(d) + ip 2 {d) + ^3(0), each item ip^{6) = J2T=-oo ^3a+n exp[i(3s + /x)0] being an 
eigenfunction of the rotation operator U r = exp[— (27ri/3)m 2 ]. Arranging the items in the form 
of a 3-component spinor (compare with eq. (|I~4"D 




*(0) = MO) (46) 
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the rotation matrix acquires the form (we set again p = 1 

'p* 
(3* 
1 




U r =(0 F |=exp<-^|0 1 0] \^exp(-^-^=\ 8 ). (47) 



On the last step we dropped a trivial phase factor. The diagonal matrix As is one of the 
standard generators of the group SU(3) and (3 = exp(27ri/3) = -1/2 + i v / 3/2. Now, since the 
factor e changes the index \i by ±1, we have in such a representation 




,±tf _ « ±ie \ ±; A+ = At = 1 ; A+A_ = /. (4* 



The matrix A + shifts each element in the column ([46p by one position down and puts the 
lowest component on the very top when the matrix A_ is the reciprocal transformation. Both 
the matrices are traceless. Therefore, the kick operator gets the form 



U fc = exp 



k 

-2 { 



-^(e ie A+ + e^A_ 



rx ] >{- l -[v(X + + X^-iv'(X + -X^]}. (19) 



In terms of the standard SU (3) generators the matrix A + reads 



A+ = i(Ax + A 4 + A 6 ) - l -{\ 2 - A 5 + A 7 ) . (50) 



The commuting matrices A± are simultaneously diagonalized, 



7? 0\ ((5* N 

x (diag) = | Q Q A (dia S ) = Q p Q \ ( 51 ) 



1/ V 1 



with the unitary transformation 




(52) 



Correspondingly, the diagonal form of the kick is 

Uf a9) = diag{e- lv+ ,e- iv -,e- lV0 ); v± = k cos(# ± 2tt/3); v = v = k cos6. (53) 

Obviously, v+ + v_ + v = 0. 

The resonant Floquet operator is now represented as 

= (54) 

where 

Ug s) = $iU r $ fe uf a9) = /3* [/ - (1 - (3)1 ■ 1 T ] Uf a9) (55) 
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with the one-column matrix I being equal to 




(56) 



Due to the factorized form (|55|) of the matrix IJ^ S \ it can be easily diagonalized. The eigen- 
vectors prove to be equal to 



^ = a - P) 



u, 



-1/2 



rrr(diag) _ 



I; n, v = +, -,0 



(57) 



while the eigenvalues u ll = e ie ^ 9 > satisfy the cubic equation 
1 - (1 - /3) I 



0. 



{hi 



After this equation has been solved, the resonant Floquet transformation is represented as 

U^ s) = exp (-z$e$ f ) = exp {-iw{6) n(0) ■ A] ; e = dia#(e + , e_, e ) (59) 
in terms of the unitary matrix $ = $ where the tilded matrix 



consists of the eigenvectors (0). On the last stage we expressed the 3x3 hermitian matrix in 
the exponent in eq.([59"D in terms of the SU(3) generators A . The coefficients are given by 



wn a = -tr [A a $e$ f ] ; a = 1,2, 



(61) 



The vector n is a unit vector in the 8-dimensional adjoint space. The transformations described 
give complete solution of the problem for q = 3 resonance. However, the final expressions in 
terms of the roots of the cubic equation ( |58| ) turn out to be too cumbersome and we do not cite 
them here. In the similar way, analytical solution can be found also for the resonance q = 8 
which is linked to the SU(4) group. However, algebraic problems increase very rapidly with q. 

Qualitatively, the situation near the lowest resonance does not essentially differ from that 
near the resonance q — 4. Figs. 8a,b demonstrate the transition from the regular motion to 
diffusion near the resonances q = 3 and q — 4. In both cases the resonant growth is seen at 
the first stage, which then changes by a sharp turnover. Further evolution crucially depends 
on the detuning. If it is below some critical value, the asymptotic plateau begins immediately 
after the turnover. However, above this value the plateau is preceded by the stage of diffusion. 
It is clearly seen that the slope at the latter stage (k • E oc t) is twice as smaller as that at the 
resonant stage (k ■ E oc t 2 ). 



5 General Case 
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5.1 Floquet Operator for a Resonance of Arbitrary Order 

The consideration presented above can be easily extended to a resonance of an arbitrary order 
q. The wave function is expressed in the form of a g-component complex vector 

y T (6) = (Md)M8)-U0))- (62) 

The component 

oo 

ijf,= ^2 C sq+f ,exp[i(sq + ^)6}; fi = 1, 2, 3, q (63) 

s=— oo 

is an eigenfunction of the operator Ur = exp(— 2itim 2 /q), which belongs to the eigenvalue 
exp [(— 2m/q) mod(fi 2 , q)]. Therefore, the rotation for the resonance (p, q) is implemented by 
the diagonal matrix 

U r (p,g)=exp(-^A ) (64) 

where all matrix elements of the diagonal matrix Ao are integers from the interval [l,q\. The 
kick operator looks as in eq. (^) where now 

A+ = A - = (°/ T o) (65) 

and is the (q — l)-dimensional zero column vector when / is the unit matrix of the same 
dimension. Both these matrices are traceless and their properties are similar to those in the 
case q = 3 described above. In particular, they are simultaneously diagonalized with the matrix 



of eigenvectors which are the discrete plane waves inside a sample of the length q in the angular 
momentum space. These waves satisfy the periodic boundary conditions at the ends of the 
sample. The diagonal representation of the kick operator is a natural extension of eq. (|53|) 

Uf a9) = diag (e- iv \e~ lv \ ...,e~ lv «) ; V/t = kcos (o + —//J ; £^ = 0. (67) 
Thus, matrix elements of the kick operator are given by the finite sums 



( u ^ = ^E e " ¥(M "^ e " rae ( 6 ^ 

while the resonant Floquet matrix reads 



(V ( ;f)^ = e-^ modi » 2 ' q) (lJ k )^. (69) 

Additional phases in eq. (|6~9D , which become quasi-random when the integers p and q are large 
enough, spoil the uniformity along the sample and disrupt the plane waves (|66|). This results 



15 



in localization of eigenvectors of the resonant Floquet matrix |5], IB]. On the other hand, 
the rotation operator is converted with the reciprocal transformation into the matrix 



7m- 



1 q 



e « e 9 . 



(70) 



Eqs. (|67| , |70|) bridges our representation to that used in the ref. |12| . 

However, the both representations yet considered do not permit extending to the case of 
a finite detuning from the resonance. To make this possible, we have at first to represent the 
resonant Floquet operator in the from of an united SU(q) transformation (compare with (|59|) ), 



Trj(res) 

p,g 



A 



cxp 



= e i exp [— i{v Xr + v 'Xj) 
exp(-i$e$ t ) = exp (—iw n A) = exp (— iR.( res n 



where 



A fl = ^(A + + A_) ,A / = i(A + -A_) 



(71) 



(72) 



when the matrices A a ; a = l,2,...,g 2 — 1 are the standard generators of the fundamental q- 
dimensional representation of the group SU (q). As before, the irrelevant phase factor generated 
by the trace trXo is omitted in final expression. This implies the condition tre = X^^m = 0- 

Diagonalization of the matrix (|69|) cannot be carried out analytically if the dimension q 
of the fundamental representation exceeds 4. Nevertheless, some generic conclusions can be 
drawn from eq. ([71]) even without explicit knowing the functions w(8) and n(0). The diagonal 
eigenvalue (quasienergy) matrix e is connected to the them as 



w 



where $ is the unitary matrix 



(73) 



(74) 



of the normalized eigenvectors of the resonant Floquet operator. Matrix elements of the q 2 — 1 
matrices $^A a $ form in (q 2 — l)-dimensional adjoint space q 2 vectors Y^'^ with the components 

0W f A a 0M = y/2T(* v \ (75) 



The diagonal part of ([73|) reads in these terms 

e M = V2w (n ■ T^) , = 

whereas the off-diagonal part yields the orthogonality condition 

n • T M = 0; [i^v. 



(76) 



(77) 
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Because of hermiticity of the generators A a , all q vectors are real. For the same reason 
the remaining vectors with // ^ v make up q(q — l)/2 mutually complex conjugate pairs Y^' u \ 
Y(">w = y^'^ . Using the known property 



1 



a \ y / 

of generators of the SU(q) group, we easily find for scalar products of the vectors T 

r M * ■ = - -^j^ . (79) 

In particular, any given vector Y^'^ with fi ^ v is orthogonal to all others including the 
vectors Y^. Therefore, the set {Y^'^; /i ^ i/} constitutes an orthonormalized basis in the 
q(q — l)-dimensional subspace of the adjoint space. As to the vectors which lie in the 
complementary orthogonal (q — 1)- dimensional subspace, they are neither orthogonal, 



Y^) . y^') 



(see (^)), nor linearly independent because of relation 

T M = tr A = . 



10) 



51) 



Excluding with the help of this relation one of the vectors, we can use the rest of them to con- 
struct an orthonormalized basis in the complementary subspace as well. In fact, the overfullness 
mentioned becomes insignificant in the limit of large q. 
On the other hand, it follows from (|73f) that 



52) 



wn = -tr (e $ f A$) = -^= ^ e^T M . 
Using also eq. (|80| ) together with the fact that ^ = we obtain finally 



w 



53) 



Note that the vector n entirely belongs to the (q — 1) -dimensional subspace spanned by the 
vectors Y^. 



5.2 Resonant Evolution 



In general case we have similar to eq. (|19|) 



AM(t) 



-<&\ —t + ie 
dv 



ic i 



$ T — e 



ict 



AM(i) • A = [M {0) t + M (1) (t)] ■ A . 



54) 
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The equation obtained is nothing but the matrix form of the result of the paper [15]. 
One can advance as follows. First of all we find from this equation 



M (°) = -itr(^$tA$ 
2 \ do 



— V f' tW 

V2^ M 



J5) 



and 



-tr e 



d0 ' 



-ict 



$t,\<M<J) 



d0 



56) 



Terms with [i — v drop out the latter sum. Therefore, the two vectors M^°) and lie in the 
mutually orthogonal subspaces. The vector M^'(i) is a quasi-periodic function of time for any 
fixed value of the angle 9. On the other hand, direct calculation gives at the moment t — 1 



AM(1) = Uj. 



-u'Aii + uAj = u(t>, f ') • A . 



(87) 



Here we took into account that the free rotation operator Ur commutes with the angular 
momentum M and that the kick operator reads = exp [— i{v Xr + v'Xj)] (see eq. (|7T|)). 
Components of the real vector u(v,v') are easily calculated as 

u a = -t>tr(A a A/) - -v'tr(\ a \ R ) . (88) 



Comparison with eqs. (84)-(|86|) taken at the same moment t = 1 gives 

d(f> {u) 



d0 



U . T (^) 

v 7 ^ ... . , — 



59) 



so we come to 



M< ' = ^(u-T(">) Y 



(90) 



and 



(91) 



instead of eqs. (|85[) and (p6|) . All derivatives disappear from these expressions. 

There exists another representation of the evolution which is more suitable for our further 
purpose. It comes from the connection 



AM(t) 



dd 
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= -w' (n • A) t - w n' T / dr e -^ (resV A . (92) 

Here 

3 res) = wn-A (93) 

and the matrix relation 

exp (iw n • A r) (a • A) exp (— iw n • A r) = a T exp (— iw n • A r) A , (94) 

where a is a (q 2 — l)-dimensional vector, has been used. The matrices A a whose matrix elements 
look as 

(A a ) 6c = -2if ahC} (95) 

with the quantities f a b c being the structure constants of SU (q) group, are generators of the 
adjoint representation of the group. From dynamical point of view, the interrelation of the 
operators 7i^ res ^ = w(n ■ A) and £( res ) = w(n ■ A) is that between Hamilton and Liouville 
operators. 

Now, let the vector x^ be an eigenvector which belongs to the eigenvalue l a of the Liouville 
matrix £( res ). In accordance with the meaning of the adjoint representation, this vector obeys 
the condition 

[n-A,x (a) -A]_ = [ a x (a) -A. (96) 

As one can easily convince oneself by direct substitution, the vectors Y^'^ satisfy the equation 
of the form (|9~6"D with the eigenvalues e„ — e M . This elucidates the meaning of the vectors T as 
the eigenmodes of the Liouville operator. It is convenient to re-number q(q — l)/2 solutions 
with /i > v with the help of the superscript a = (3 = 1, 2, q(q — l)/2 so that 

T<*"> e u -e^l p . (97) 

Then for the solutions with v > fi 

X { -? ] =X {0) *i U = -h- (98) 

There also exist q — 1 real eigenvectors a = 1,2, ...,q — 1 with zero eigenvalues l a = 0. All 
these zero modes are pairwize orthogonal linear superpositions of the vectors Y^. Obviously, 
one of them is just the vector = n defined in eq. ([83]). 

Again, we can separate in eq. (^) the linearly growing and quasi-periodic contributions 
and then exclude the derivatives w' and n' by comparing the result with eq.(^). In such a way 
we arrive at 

M<°> = ^ (u • X {a) ) X (a) , (99) 

a 

M(1, ( ^ > = £!^( u 'X (e, )x (e| •. ( 10 °) 

The result fllOOD is just identical to ©. As to eqs. © and ([lOOD , they express the vector 
M 1 - ) in terms of different sets of vectors which are linear combinations of each other and belong 
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to the same set of the degenerate zero eigenvalues. Thus we come to the conclusion that the 
resonant growth entirely originates from the zero Liouville modes while the non-zero ones yield 
quasi-periodic evolution. 

Kinetic energy at the moment t is equal to 

E(t) = ^(AM a (i)AM t (t))$jA a A 6 * = ^1a([AM a (t)f). (101) 

a, b a 

The symbol (...) stands for the ^-averaging. As before, the initial state \l/o is chosen to be 
isotropic, \Ed = (0, 0, 1), so that ^IXaX^o = S a bf] a with 

Va = if a < (q - l) 2 - 1; rj a = 1 if (q - l) 2 < a < q 2 - 2; rj^ = 2(1 - l/q). (102) 

The first set of indices enumerates such matrices A a which include the SU(q — 1) generators. 
Such matrices annul the initial state \l/o- The second set marks 2q — 1 non-diagonal A-matrices 
with non-zero elements in the qth columns and rows. The matrix X q is the last diagonal 
generator of the group SU(q). So, only those modes are significant which have appreciable 
projections onto the (2q — l)-dimensional active subspace indicated above. 
According to eq. ( |101| ), the energy resonant growth rate is is given by 

r(p } q;k)= 1 -Y,Va([M^} 2 }. (103) 

a 

Contribution of the non-zero modes fluctuates with time and approaches asymptotically a 
finite positive value 

Xoo (p, q; k) = hm l - £ Va ( [M« (t)] 2 ) . (104) 

a 

Estimation of the interference terms is less certain. If the stationary points dlp/d9 = are not 
degenerate, terms oc \ft arise together with the linearly growing ones. 

5.3 Vicinity of a Resonance 

Motion near the point of a resonance of the order q is described by the quasi-Hamiltonian 

H = KH (res) + k 2 Q(k) = KH {res) + k 2 Q (0) + k 3 Q (1) + ... 



F o (0) + - {J, Fi(0)} + + -JF 2 (0)J + ... (105) 



where 

1 = 1® J, ¥j(9) = Fj{9) + Fj(9) ■ A . (106) 

The operators Fj are found from the matrix analogs of the conditions (^), (|T0|) . For the first 
correction eq. @ gives F 2 (9) = I while 



I drF 1 (9;-r) = - [ drAI(-r) 
Jo Jo 



(107) 
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f dr¥ (8;-r) = - f dr ([AJ(-r)] 2 + {AJ(-r), Fx(0; -r)} + ) . 
Jo Jo 



'108) 



Here 



Wj(9;-t) = e- tnv,T ¥ j {9)e %nv ' T , AJ(-r) = «AM(-r). (109) 
The l.h.s. in eqs. ( |107| ), ( [L08| ) are calculated with the help of the connection 



[ dT¥j(8; —t) = Fj(8) + Fj(8) fdre* e " T • A. (110) 
Jo Jo 

Because of zero modes, a reciprocal to the operator £( res ) does not exists. This prohibits 
from straightforward integration over r. To avoid the said difficulty, we regularize operator 
^(res) _ £( res ) _ w ] iere infinitesimal 5 must be set to zero at the very end of calculation. 
Then the operator 



1 1 iC ( - rea> 



dre^i v = — __ (HI) 



becomes well defined. Now, in virtue of eqs. fl92|) , (|99|) and (|100|) 



AJ(-r) = k u T 1 ~ e% S I ■ X . (112) 

1 - e-< 



Substituting eqs. ( P-ldj ), ( P-l 1| ) and ( |112j ) in the condition ( p.07| ), we finally find 



1 + - e i£ - { r s) 

'■"- ^JA.-,-^ ^ (113) 



Due to skew symmetry of the Liouville matrix, (£( res )) T = _ £( res ) ; the matrix Fi(#) is hermitian 
as eq. (|105D implies. Other conditions are solved in the similar way. However, corresponding 
expressions are very bulky in general case. When the kick parameter k ^> 1 the typical number 
of harmonics in operators Fj(0) is proportional to k, with a rapidly increasing coefficient £(?). 
This results in very fast diminishing of the widths of resonances when q grows. 

In the approximation (|105|) the quasi-Hamiltonian TC is formally equivalent to the Hamil- 
tonian of a quantum particle with q intrinsic degrees of freedom which moves in a (q 2 — 1)- 
component inhomogeneous "magnetic" field. In a sense, such a motion is a quantum analog of 
the classical phase oscillations near a non-linear resonance. 



6 Convergency Problem 

Consideration presented above shows that in some domain of the detuning from a quantum 
resonance of a finite order q evolution of QKR looks like a conservative motion described by 
the effective time-independent quasi-Hamiltonian with a discrete quasi-energy spectrum. A 
few lowest terms of the expansion ( |105| ) allow to predict with great accuracy the evolution 



for a very long time. More than that, in the range An of such a domain (the width of the 
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resonance), which is determined by the condition that the influence of higher corrections is 
week, the accuracy is improved with each further correction kept. Nevertheless, it is clear that 
the formal expansion ( |105| ) cannot converge. Indeed, within the width of any strong resonance 
of a relatively small order there exist an infinite number of resonances of large orders which give 
rise to unrestricted resonant energy growth. Independently of the number of corrections taken 
into account, the quasi-Hamiltonian approach fails to reproduce such a growth which implies 
continuous spectrum. 

Actually, the resonant rate r(p, q; k) decreases when the order q increases and becomes 
exponentially small if q noticeably exceeds the typical localization length. In Fig. 9 we plot the 
empirical dependence of the resonant growth rate on the resonance order q. At each q the ratio 
p/q is chosen to be as close to the "most irrational" number (v^5 — l)/2 as possible. Our data 



are in agreement with earlier results from |12|, [18], [T{|. The solid line shows the fit with the 
help of the semi-empirical formula 



2k 2 

r(p, q;k) = — exp {-q/2l) (114) 
3q 



proposed in pi 



The qualitative arguments presented in p0| connect the exponential suppression of the 



resonant rate with the localization and tunneling in the momentum space. The wave function 
corresponding to an eigenvector (f>^ reads [|i~2] 



oo q 1 s=oo 

^\9)= ^ ] e im(> = Y^\0o)e lv{e - eo) — e^ (e ^ o) . (115) 



2tt 



One must distinguish here between the coordinate eigenvalue 6q an d the argument 9 of the coor- 
dinate representation. In the angular momentum representation this equation yields IV^+sJ = 

|0^(^o)|> so that the angular momentum distribution is periodic with the period q. 

On account of the quasi-random phases in the matrix ( |S9"D (see also discussion below this 
formula) the g-dimensional eigenvectors (f>^' are, in fact, localized so that the number I of 
their appreciable components is much smaller than q. Then overlap of neighboring bumps of 
an eigenfunction ipm^ is, typically, exponentially weak. This resembles the eigenfunctions of 
a particle moving in a periodic chain of potential wells. The resonant grows of the angular 
momentum of QKR is similar to the transport through the chain of the particle which initially 
was localized in some well. The rate r(p, q; k) is an analog of the mean value of the squared group 
velocity of such a particle. The latter is proportional to the exponentially small probability of 
tunneling between neighboring wells. This interpretation is in agreement with relation |12| 



r(p,q;k)= 1 -J2^) 2 l^l 2 ) ( 116 ) 



which directly follows from eq. (pi ) and contains the weighted-mean value of the squared 
"group velocities" e' . Returning to the expression ( |103| ), we conclude that in the case q ^> / 
the Liouville zero modes with exponential accuracy lie in a subspace orthogonal to the active 
one. 

The aforecited arguments show that in the case of large q the resonant growth reveals itself 
only on a very remote time asymptotics owing to the tunneling between localized parts of 
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the globally delocalized quasienergy eigenfunctions. In particular, if such a resonance hits the 
domain of a strong one, it happens after the exponentially large time 



A 



nr(p, q; k) 



1/2 

ocexp(g/4Z) (117) 



when quadratically growing contribution becomes comparable with the height of the plateau 
E p i p» A/k formed due to the influence of the strong resonance (see Figs 7a,b). We found it 
difficult to calculate the coefficient A analytically but numerical simulations show that under 
condition k 3> 1 it weakly depends on the kick parameter k as well as on the order of the strong 
resonance. 

Exponential effects of such a kind, which are characteristic of the tunneling, are well known 
to be beyond the reach of perturbation expansion. For this reason they cannot be described in 
the framework of the quasi-Hamiltonian method. This approach reproduces only those features 
of the motion which are determined by the discrete component of the quasienergy spectrum. In 
particular, contribution of the non-zero modes much faster attains its asymptotic value (|104| ). 
As a result, for all weak resonances inside the width of a strong one the time dependence of of 
the part MW(t) is dictated during exponentially long times t < t r by their strongest brother . 

On the other hand, if the order of the resonance we are interested in is very large, q 3> I, 
and this resonance lies in the region of typical irrationals being far from all strong resonances, 
already a very small detuning suffices for killing the quadratic growth with exponentially small 
rate. At the same time, such a shift does not influence the term MW(t) which reproduces 
on exponentially large (though finite) time scale all characteristic features of the "localized 
quantum chaos" [pi]]. 



The behaviour is most complicated and ambivalent in the transient region q ~ I. In this 
case a number of resonances with comparable and moderate orders are neighboring and their 
domains can overlap. The expansion near one of them forms a plateau which lasts until the 
quadratic growth in a next resonance of the same strength reveals itself so that the original 
expansion fails. However, the expansion near the new resonance cuts off the growth and forms 
a higher plateau until a next resonance comes to the action. Such a pattern of repeatedly 
reappearing regimes of resonant growth has been discovered in [0 . 



7 Summary 

In this paper we propose on the example of the QKR model the concept of the time independent 
quasi-Hamiltonian of a quantum map. The regimes of quantum resonances, which take place 
under conditions k = const and <; = T/4tc = p/q, play a crucial role in our construction. 
The motion in the very point of a quantum resonance with the order q is, generally, exactly 
described by a continuous transformation from the SU (q) group. The two qualitatively different 
contributions: growing with time and saturating, in the evolution of the QKR come respectively 
from the zero and non-zero modes of the generator of the SU(q) transformation in the adjoint 
representation of the group. A perturbation expansion exists near the point of a given quantum 
resonance, which provides quite a good description of the motion within some domain - the 
width of the resonance. Inside the width of a strong resonance the motion is mastered by the 
resonance. This motion is proved to be similar to that of a quantum particle with q intrinsic 
degrees of freedom along a circle in an inhomogeneous (q 2 — l)-component "magnetic" field. 
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The width of a quantum resonance strongly depends on its order q. The resonances with 
smallest orders are the strongest ones and have maximal widths. The motion within the width 
of such a resonance, being dominated by it, is proved to be regular. In all cases save the 
two boundary resonances q = 1,2 the regular quantum motion exists in spite of the fact that 
the corresponding classical motion is chaotic and exponentially unstable. The widths of these 
resonances vanish in the classical limit k — > oo , T — > , K c = kT = const. Such a situation 
holds as long as the condition q <C / takes place, where / is the localization length for close 
typical irrational q, although the widths rapidly diminishe with q. In the opposite case of very 
large orders q ^> I and p ~ q the motion weakly depends on the q and as well as on the detuning 
k. The motions with rational and irrational q differ only on a very remote time asymptotics. 
During a long though finite time the motion reveals universal features characteristic for the 
localized quantum chaos. 
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Figure 1: The kinetic energy < J 2 /2 > versus 
the number of kicks t for the main resonance 
q — 1. The dotted and solid lines show pre- 
dictions of the proposed theory in the first and 
second orders, respectively. Open circles cor- 
respond to the exact QKR map. 
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Figure 3: The kinetic energy < J 2 / 2 > versus 
the number of even kicks for q = 2. The solid 
line shows evolution predicted by the quasi- 
Hamiltonian (|28|) while the open circles cor- 
respond to exact numerical simulations. The 
dotted curve at the bottom indicates the devi- 
ation of the theory from the exact solution. 
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Figure 2: The kinetic energy < J 2 /2 > ver- 
sus the number of kicks t for q = 2. Dots 
and crosses show predictions of the theory in 
the first and second orders. Open circles cor- 
respond to the exact QKR map. The two (al- 
most) symmetric branches are due to the only 
even or only odd kicks respectively. 
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Figure 4: The kinetic energy < J 2 /2 > versus 
the number of even kicks for q = 2. The solid 



line results from the classical limit (|30D ; crosses 
and open circles correspond to the exact QKR 
map at the same effective classical parameter 



K e but different k. 



26 



1 1 1 1 1 1 1 


o t=4n 


- k=10 


• t=4n+l 




» t=4n+2 


*** 


- t=4n+3 










r / 




X(oo,k=10) 









40 80 120 160 200 

t 

Figure 5: The four branches of the function 
x(t, k). Asymptotically all of them converge to 
the same limit k) shown by dashed line. 
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Figure 7: The plateau height E p l versus de- 
tuning k near the resonance q — 4. The 
solid line corresponds to the theoretical rela- 
tion Epl oc l//c. 
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Figure 6: The kinetic energy E versus the kick 
number t for q — 4. The solid line shows 
the second order perturbation theory and the 
points correspond to the exact QKR map (each 
4th kick is kept in the main part and each 
500th in the inset). Deviations of the theory 
from the exact map is indicated at the bottom. 
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Figure 8: The plateau height E p l versus de- 
tuning k near the resonance q — 3. The 
solid line corresponds to the theoretical rela- 
tion Epl oc 1/k. 
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Figure 9: The crossover region near the reso- 
nances q = 3 and q — 4. The upper curves 
correspond to k outside the resonance widths. 
The diffusion is clearly seen after a short ini- 
tial resonant stage while at smaller values of 
the detuning the diffusion does not develop. 
Double log scale is chosen to show all stages of 
the time evolution. 




400 600 

q 



1000 



Figure 10: The resonant growth rate r(q) ver- 
sus the resonance order q. At each q the p/q 
ratio is chosen to be the closest to the most 
irrational number (\/(5) — l)/2. 
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